Limits on inelastic dark matter from ZEPLIN-III 
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We present limits on the WIMP-nucleon cross section for inelastic dark matter derived from the 
2008 run of ZEPLIN-III. Cuts, notably on scintillation pulse shape and scintillation-to-ionisation 
ratio, give a net exposure of 63 kg • days in the range 20-80 keV nuclear recoil energy, in which 6 
events are observed. Upper limits on signal rate are derived from the maximum empty patch in the 
data. Under standard halo assumptions a small region of parameter space consistent, at 99% CL, 
with causing the 1.17 ton • year DAMA modulation signal is allowed at 90% CL: it is in the mass 
range 45-60 GeV with a minimum CL of 88%, again derived from the maximum patch. This is 
the tightest constraint on that explanation of the DAMA result yet presented using a xenon target. 



PACS numbers: 95.35.+d, 29.40.Mc, 29.40.Gx 

Dark matter in the form of weakly interacting massive 
particles (WIMPs) which scatter predominantly into a 
higher-mass state has been proposed [ij as an explana- 
tion of the annually modulated event rate in DAMA/Nal 
and DAM A/LIBRA 01 which is also consistent with the 
upper limits on WIMP-nucleon elastic scattering rates 
from other experiments 043 ■ In such inelastic dark mat- 
ter (iDM) models, scattering with energy transfer due 
to a WIMP of ground state mass and mass change 5 
requires a minimum relative speed 

V2ton£'r ( Mn ^ ) ' ^ ^ 

where ttin is the nucleus mass and /iN is the reduced 
mass of the WIMP-nucleus system. A non-zero 6 results 
in a recoil spectrum that is zero at low energy and more 
sensitive, compared with elastic scattering, to the upper 
tail of the WIMP velocity distribution. WIMPs with 
velocity below (2(5//in)°^ will not scatter inelastically at 
all and so, for a given local escape velocity, more m^- 
d parameter space is accessible to heavier target nuclei. 
However, systematic uncertainty in the expected relative 
rates in different targets due to nuclear form factors and 
WIMP velocity distributions grows with the difference in 
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atomic mass [6]. On balance, xenon is well suited to test 
iDM models that would, by predicting a modulated rate 
of scattering against iodine nuclei, explain the DAMA 
observation. 

ZEPLIN-III (described in detail in Refs. 0, S]) is a 
liquid/gas detector designed to search for WIMPs scat- 
tering against xenon nuclei in the 6.5 kg fiducial liquid 
volume. It is built of low-radionuclide components, en- 
cased in hydrocarbon and lead shielding, and operated in 
the Palmer Laboratory at Boulby Mine beneath 2850 m 
water-equivalent rock overburden. Events are charac- 
terised by two light signals recorded by an array of 31 
photomultiplier tubes (PMTs) . The summed scintillation 
signal from the liquid is denoted by 5'1. A 3.9 kVcm~^ 
electric field in the liquid extracts ionisation charge from 
the interaction site, drifts it to the surface and forces 
emission into the gas layer above; there, an electrolumi- 
nescence signal, S2, is produced. As described in Ref. 
events with one 5*1 and one 52 signal were selected and 
cuts made, based on the pattern of light distribution, to 
remove multiple-scintillation, single-ionisation events. 

An event's electron recoil equivalent energy, denoted 
by Ecc and measured in keVee, is derived from the pulse 
area of the 5*1 signal, normalised to 122 keV photoab- 
sorption using a '^^Co 7-ray source. Discrimination be- 
tween nuclear and electron recoil events is achieved pri- 
marily through the ratio of scintillation and ionisation 
signals. Additional discrimination has been achieved here 
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Figure 1: The eflFect of the n cut on the 20-25 keVee bin, over 
all S2/S1. The points and right scale correspond to AmBe 
data, fitted by a gamma distribution. The outlined histogram 
and left scale correspond to ^''^Cs data before the cut with 
a shaded region corresponding to the 13% of electron recoil 
events which are not rejected by the cut. 



using scintillation pulse shape. Recoiling electrons and 
nuclei produce different proportions of the singlet and 
triplet excited dinier states, which have lifetimes of 4 and 
22 ns respectively [£] . PMT traces are sampled at 2 ns in- 
tervals and so the mean arrival time of the 51 pulse area, 
denoted by ri, is a useful discriminator. The timing of 
AmBe neutron calibration events within each 5-keVee bin 
from 5 to 40 keVee is well described by gamma distribu- 
tions in 1/ri [lO|. Fitting a polynomial in Eee to the 
medians of the gamma distributions produces a cut on 
Ti with 50% signal acceptance. Fig. [T] shows the separa- 
tion of the two recoil types and the effect of the cut in an 
example bin. The power of the timing cut to reduce elec- 
tron recoil background increases with energy, as seen in 
Fig. [21 mainly due to a narrowing in the ti distribution 
of electron recoil events. 

AmBe calibration data were also used to obtain the 
52/51 distribution of elastic nuclear recoil events which 
pass the timing cut, as a function of Ece- As in Ref. 
the logiQ (52/51) distribution was fitted by a Gaussian in 
each energy bin, and the energy dependence of the fitted 
means and standard deviations parametrised by a power 
law to define a cut with 47.7% signal acceptance. Charge 
recombination causes 52 and 51 to be microscopically 
anticorrelated at a given energy; in principle, therefore, 
the 52/51 distribution has some dependence on the recoil 
energy spectrum. However, the low level of field-induced 
51 suppression observed for nuclear recoils in xenon [11] 



Figure 2: The efficiency of the ti cut on calibration events 
over all S2/S1, for energy bins 5-40 keVee. Red triangles are 
AmBe data and blue squares are ^^''Cs data. 



suggests that the effect is relatively small. Here we have 
assumed, as xenon experiments historically have, that the 
52/51 distribution at fixed 51 for neutron calibration 
events is an adequate approximation to that for signal 
events. After efficiencies from dead time, pulse-finding, 
event reconstruction and the cuts on 52/51 and ti, the 
net exposure for signal events is 63 kg • days, with 5% 
uncertainty due to neutron calibration statistics. 

Nuclear recoil-equivalent energy, .Er, is determined as 
in Ref. [4] from Ecc via a conversion factor: 
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where 5o and 5n are the field-induced suppression factors 
for the light yield of electron and nuclear recoils and L^s 
is the zero-field light yield of nuclear recoils relative to 
that of electron recoils. An energy range of 20-80 keV nu- 
clear recoil energy (8.4-38.3 keVee) was chosen to include 
the majority of events predicted by the quenched, inelas- 
tic WIMP-iodine scattering interpretation of the DAMA 
modulation [12]. 

Fig. [3] shows the six search events which passed all cuts. 
The combined efficiency of the cuts on 52/51 and ri for 
search data (6 in 1.3 x 10^) is no higher than for electron- 
recoil calibration data (7 in 8.5 x 10'*), suggesting that 
the surviving search events may well constitute the tail 
of the electron-recoil background population. Without 
the timing cut the box would have contained 27 events. 

For WIMPs which couple equally to protons and neu- 
trons, the differential rate for spin-independent WIMP- 
nucleus scattering in a target of total mass Afx is given 



Figure 3: Search data passing all cuts except those on Ecc and 
S2/S1. Events passing all cuts are highlighted by triangles. 
The solid and dashed lines show the mean plus and minus two 
standard deviations for elastic AmBe calibration events, and 
the vertical lines indicate 20-80 keV nuclear recoil equivalent 
energy. 



by: 



dR 



[E^,t) = ^I^A^F\q)r SJ-^, (3) 



where is the local WIMP density, A is the atomic 
number of the target nucleus, an is the WIMP-nucleon 
cross section, /!„ is the WIMP-nucleon reduced mass and 
/ (w, t) is the WIMP velocity distribution in the target 
frame. A Helm form factor was used: 
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for momentum transfer q, where the effective nuclear ra- 
dius is taken to be r.r, AAA 3 — 5 fm, the skin depth 
s = 1 fm and ji is a spherical Bessel function. 

Recoil energy spectra were calculated under a standard 
halo model: = 0.3 GeVc~^ cm"'^, a Maxwellian veloc- 
ity distribution with vq — 220 ms^"'^ truncated at escape 
velocity Vosc in the galactic frame, and an Earth velocity 
parametrised as in Ref. The underlying spectrum for 
given , S and (T„ was modified by the energy resolution 
and efficiency of ZEPLIN-III and then averaged over the 
83-day run to produce a signal model. The energy reso- 
lution, dominated by Poisson statistics of photoelectron 
production and the variance of the single-photoelectron 
response, is ct/E'r = 1.5 (jj^R./keV) 

The maximum patch statistic |l4| was used to derive 
single-sided upper limits on the rate of signal events in 



Figure 4: (color online) 90% limit on cr„/cm^ as a funciton of 
mass and splitting. The upper left region predicts no inelastic 
scattering during the run. 



the 20-80 keV range. No background estimate is used; 
consequently, the null hypothesis cannot be ruled out by 
this method. Events were mapped onto a plane of uni- 
form signal density by integrating the signal spectrum in 
-Er and the fitted profile in 52/51. For models in the 
previously un-excluded region of iDM parameter space, 
the largest empty rectangle in the re-mapped search box 
has a fractional acceptance of 0.73-0.75; this implies a 
90% CL limit of 5.4-5.1 expected signal events in the 
box. The resultant limits on (t„ for Vcsc = 550 kms~^ 
are plotted in Fig. 01 

Signal modulation spectra for the combined DAMA ex- 
periments were constructed with resolution as described 
in Refs. [l^[lB| and parametrised in Ref. [13 ■ An iodine 
quenching factor of 0.08 [3 113 used; the exclusion 
results are relatively insensitive to channelling effects |20| 
which are, conservatively, omitted. The parameters m^, 
5 and cr„ were fitted, by minimizing to the observed 
modulation amplitude in 0.5-keVee bins from 2-10 keVee 
and a single 10-20 keVee bin, following Ref. [13 ■ A 90% 
confidence interval for the local escape velocity from Ref. 
[21] is 498-608 km s~^ and the cross section excluded by 
ZEPLIN-III depends on the true Vcsc- Non-Maxwellian 
velocity distributions would cause a similar systematic ef- 
fect. Fig.[5]shows the ZEPLIN-III constraints on parame- 
ter space consistent with causing the DAMA modulation 
for three values of Vgsc- DAMA-explaining cross sections 
are excluded at the 88% confidence level. Fluctuations 
of ± 1 • CT in the cut efficiencies derived from neutron cali- 
bration would change this minimum CL within the range 
87-90%. 
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Figure 5: In m^-S space, the confidence level at which ZEPLIN-III excludes the lowest value of (t„ consistent, at 99% CL, with 
causing the DAMA modulation. Three values of Vcsc are shown: (from left) 500, 550 and 600 km s~^. 



In summary, a search of 63 kg • days net exposure with 
a xenon target yielded 6 candidate events in the range 
20-80 keV nuclear recoil equivalent energy. They were 
consistent, both in number and scintillation-to-ionisation 
ratio, with belonging to the tail of an electron recoil 
background population. Single-sided upper limits were 
set on the WIMP-nucleon cross section, constraining the 
DAMA-explaining region of iDM parameter space: for a 
standard halo model there remains a 90% CL allowed re- 
gion for WIMP masses in the range 45-60 GeV c~^, with 
minimum CL 88%. This is more stringent than limits 
from other xenon and germanium experiments [^,[22, 23 1 
and supports previous exclusions fl7'| based on CRESST- 
II data. In particular, a target element of similar mass to 
iodine reduces systematic uncertainty due to the WIMP 
velocity distribution. 
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